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The layered lithium borocarbide LiBC, isovalent with and structurally similar to the supercon- 
ductor MgB2, is an insulator due to the modulation within the hexagonal layers (EC vs. B2). 
We show that hole-doping of LiBC results in Fermi surfaces of B-C pa character that couple very 
strongly to B-C bond stretching modes, precisely the features that lead to superconductivity at Tc — 
40 K in MgB2. Comparison of Lio.sBC with MgB2 indicates the former to be a prime candidate for 
electron-phonon coupled superconductivity at substantially higher temperature than in MgB2. 



The discovery of '~40 K superconductivity in MgB2 ||] 
has led to a resurgence of interest in layered compounds 
that are metallic, or may be made metallic by doping. 
MgB2 itself is, to date, a unique system, one that a num- 
ber of studies of the electronic structure confirm 
contains strongly covalent bonds similar to graphite yet is 
a good metal with strong electron-phonon (EP) coupling. 
A variety of evaluations of the EP coupling indicates that 
it is strong enough to account for superconductivity in 
the 40 K range, [pp[pT|-p^ and the observed B isotope 
shift of Tc [Q confirms that pairing is due to phonons. 
So far, attempts to make related, nearly isoelectronic di- 
borides with the same structure have been fruitless. Be 
"diboride" (actually BeB 2. 75) forms a complex compound 
[ p^ distinct from the MgB2 type structure that has been 
presumed for several band structure studies. |^,^,^ Ca, 
Sr, and Ba have resisted incorporation into the MgB2 
structure. Al can be alloyed in, resulting in a decrease 
in Tc and a subsequent collapse of superconductivity al- 
together. Alloying with monovalent ions for Mg should 
increase Tc (based on rigid band concepts) but attempts 
to do this also have been unsuccessful. 

A closely related compound that has not received at- 
tention since the discovery of MgB2 is LiBC. Although 
apparently straightforward to make , [|l6| it has been stud- 
ied very little to date. Its structure [[16|-|18| can be derived 
from the "fully intercalated graphite" structure of MgB2 
by Mg Li, and B2 BC, with the hexagonal BC lay- 
ers alternating so that B is nearest neighbor to C along 
the c axis as well as within the layer. Like the compound 
CuBC suggested by Mehl et al. it is isovalent with 
MgB2: Li has one less valence electron than Mg, but C 
has one more than B. The main features of the electronic 
structures, which are dominated by in-plane B-B and B-C 
bonding respectively, should be expected to have strong 
similarities. LiBC is insulating, however, with a gap that 
(we show below) arises from the replacement B2 — » BC, 
and not from the alternating stacking. An important fea- 
ture of Lia;BC is that the Li content can be varied within 
the same crystal structure: while this system forms stoi- 
chiometrically (a;=l), the Li content can be reduced and 
a structural change was seen only for a;=0.24. 

Ramirez et al. have presented semi-empirical 

Hartree-Fock (INDO) results for a (BC)~ layer repre- 



sentiirg LiBC in the absence of interlayer coupling. They 
obtained a band gap of 4.3 eV (acknowledged to be an 
overestimate due to neglect of correlation) , with its min- 
imum value at the K point of the Brillouin zone between 
bonding and antibonding tt [pz] bands. B and C charac- 
ter of the bands were not distinguished in their results. 
The cr bands {sp^Py hybrids) however lie above the pi: 
valence bands in a region around the F point. Albert in- 
cluded LiBC in a discussion of Li borocarbides, most 
of which are structurally very different from LiBC. 

Lr this paper we present results for the electronic struc- 
ture and EP coupling of Lia;BC (x — 0.5, 0.75, 1.0) and 
compare it with that of MgB2. We determine that once 
it is hole-doped, the states at the Fermi level include a 
bands with cylindrical a Fermi surfaces that are strongly 
coupled to the B-C bond-stretching mode (more strongly 
than in MgB2). At a; ~ 0.5, the density of states (DOS) 
N(Ei?) is comparable to that of MgB2 while the EP cou- 
pling is much stronger. These features suggest that hole- 
doped LiBC is a very favorable candidate for supercon- 
ductivity, possibly at much higher temperature than in 
MeB.,. 



FIG. 1. Crystal structure of LiBC, a direct generalization 
of that of MgB2 , showing the alternating layers of hexagonal 
BC. Li ions lie in the interstitices. 

The layered structure of LiBC, pictured in Fig. 1, has 
been determined by Worle et al. p6|] from golden hexag- 
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onal platelets (and red in transmission, hence insulat- 
ing) to be P63/mTOc (D^^, No. 194), a= 2.752 A, c = 
7.058 A= 2x 3.529 A. This structure is similar to that of 
MgB2, except that a is 9% smaller; the interlayer spacing 
is nearly identical. The BC layers are stacked in alter- 



nating fashion, so there is no B-B or C-C bonding in 
the structure. To expose the effects of A-B versus A-A 
stacking of the layers, we also consider the A-A stacking, 
which has space group P6m2 (Dg^, No. 187). 
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FIG. 2. Band structures of: (a) LiBC in the observed structure; (b) LiBC in fictitious A-A stacking structure, where 
additional dispersion produces a semimetallic band structure; (c) LiB2, the virtual crystal analog; (d) Lio.sBC in virtual 
crystal approximation; (e) C pa "fatbands; (f) B pa "fatbands; (g) C pir "fatbands; (h) B p-K "fatbands. 



We have used two electronic structure methods. A 
full-potential nonorthogonal local-orbital scheme |^l[] was 
used to obtain accurate band structures and the total en- 
ergies that we report. Li(2s, 2p, 3d), B (2s, 2p, 3d) and 
C (2s, 2p, 3d) states, respectively, were chosen as the ba- 
sis set. All lower lying states were treated as core states. 
The Li 2p and 3d states as well as the B and C 3d states 
were taken into account to increase the completeness of 
the basis set. The spatial extension of the basis orbitals, 
controlled by a confining potential |Q (r/ro)^, was op- 
timized to minimize the total energy. In addition, the 
LMTO-47 code |Q was used for the virtual crystal cal- 
culations and to obtain the deformation potentials that 
we report below. Atomic sphere radii of 3.00 a.u. (Li) 
and 1.73 a.u. (B and C) were used. 

Stoichiometric Compound. The band structure of 
LiBC (observed structure) along high symmetry lines in 
the Brillouin zone (BZ) is shown in Fig. 2(a). Valence 
bands are separated from conduction bands by a calcu- 
lated gap of 1.0 eV, with the valence band maximum 
occurring at F, and the conduction band minimum oc- 
curring at H. As in other semiconductors, this LDA value 
is likely to be an underestimate of the experimental gap. 
It will be shown below that the band structure of LiBC 
is closely related to that of MgB2, so the first question 
to address is the origin of the band gap at the Fermi 
level in LiBC, with the two possibilities being the alter- 
nating stacking of B-C layers along the c axis, and the 
alternation of B and C within the layer itself. 



This question is answered by performing two calcula- 
tions, one for a isoelectronic virtual crystal "LiS2," where 
B is an average atom between B and C with nuclear 
charge Z=5.5, and another for an alternative LiBC struc- 
ture in which the stacking of layers along the c axis is A-A 
stacking rather than the observed A-B stacking. (This is 
the structure for CuBC considered by Mehl et al. |p9[). 
These results are also shown in Fig. 2(b)-(c). The bands 
of LiB2 are isomorphic to those of MgB2, as expected, 
as it is structurally identical and chemically similar. The 
bands look different from those normally shown for MgB2 
because we have used a doubled cell (along c) to allow 
direct comparison with the bands of LiBC. The a and tt 
bands (see below) appear twice due to BZ folding, with 
the former separated by only a few tenths of eV due to 
their layered character {s — Px — Py) and the latter (p^) 
separated by 3-4 eV around and below the Fermi level 

The A-A stacking of BC layers [Fig. 2(b)] produces 
a substantial change of the electronic structure: it does 
not change the separation of the valence and conduc- 
tion bands, but it does change kz dispersion sufficiently 
that the gap disappears - it becomes a semimetal with a 
calculated band overlap of 1.7 eV. The calculated total 
energy of the A-B stacking is lower than A-A stacking (at 
constant structure, the experimental one) by 35 meV per 
formula unit, consistent with the observed A-B stacking. 

To understand the electronic structure of LiBC and 
its similarity to that of MgB2 more clearly, we display 
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in Fig. 2(e)-(h) the "fatbands" plots showing C and B 
a and tt character. The a bands are as in MgB2, with 
C character concentrated more heavily at the bottom of 
the bands (-6 to -8 eV below the gap) because the on- 
site energies Ss and Ep are lower for C than for B. The 
gap in LiBC, also evident in the DOS plots in Fig. 3, 
lies between bonding and antibonding combinations of tt 
bands that are forced by symmetry in MgB2 to connect, 
and thus enforce metallicity. The bonding tt bands have 
a preponderance of C Pz character, while the antibond- 
ing bands just above the gap are more strongly B pz in 
character. The ionic character of the B-C layer compared 
to B-B helps to push the tt character downward with re- 
spect to the tJ character, resulting in a larger proportion 
of (7 character [Ncr(e)] than in MgB2, a point of some 
importance below. 




FIG. 3. The density of states of LiBC (solid line) compared 
to that of MgB2 (dashed line, with the vertical line denoting 
the Fermi level); the gaps arise from in-plane (B-C) modu- 
lation. The main valence band peaks have been aligned to 
illustrate the similarity. 

Sub stoichiometric LixBC. In a band semiconductor 
such as LiBC, heavy hole doping will make Lia;BC metal- 
lic. If the a bands can be populated by a sufhcient density 
of holes, then by analogy with MgB2 the system should 
become superconducting, and perhaps a very good su- 
perconductor. There are at least two crucial questions. 
First, can the Fermi level be moved far enough into the 
valence bands to give a value of N(E^) large enough to 
be interesting - comparable to that of MgB2, say? Sec- 
ond, are these a states strongly coupled to the B-C bond 
stretching modes as they are in MgB2? 

We have used supercells (doubled within the layers) to 
study the x — 0.5 and x = 0.75 concentrations. For the 
X = 0.5 case, two configurations were studied to assess 
the effect of the placement of missing Li atoms. In both 
configurations half of the Li atoms were removed from 
each layer, in one case along a line parallel to the c axis, 
and in the other case in staggered positions. The differ- 
ences are very minor. In Fig. 4 the resulting DOS of 
these systems are compared with those for x — 1. The 
Fermi level moves into the valence band, as it must, and 
the behavior in the region of the gap and Fermi level is 



not far from rigid band. It is noteworthy that it is the a 
bands that are occupied first as holes are doped in [Fig. 
2(a)]. The value of N(Ef) = 0.68 eV"! at a; = 0.5 is 
nearly equal to that of MgB2. 
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FIG. 4. Densities of states of Li^jBC, for x=l, 0.75, and 0.5, 
displaced upward consecutively for clarity. The 1=0.5 curve is 
averaged over two possible supercells (see text). The main va- 
lence band peaks have been aligned (but the gaps also remain 
aligned), so the Fermi level (vertical lines) moves downward, 
roughly in rigid band fashion, as hole doping proceeds. 

Estimation of Tc- EP coupling (strength A) in MgB2, 
which drives its superconductivity, is dominated by cou- 
pling of B-B bond stretching modes to the a cylinder 
sheets of Fermi surface. [^,| 10 -1^ We calculate the de- 
formation potential of these states for the zone-center 
B-C bond stretching mode in Lio.sBC to be I? = 18.5 
eV/A, more than 40% larger than in MgB2 (where it 
is 13 eV/A 1^), resulting from the stronger B-C bond- 



ing compared to B-B. We calculate lue^ 



68 meV 



for Lio.sBC. We can obtain a realistic estimate using 
g A cx N^(Ef) jMu? and scaling from calculated 
quantities and experimental T^ for MgB2. Defining the 
ratio TliK) = Kna ^BC / KMgB2 of any property K, 
our calculations give 'R{N„) — 1.3, 7^(2?^) = 2.0, and 
TZ{ujE2g) = 1-17, finally giving 7^(A) = 2.1 as the en- 
hancement of A in Lio.sBC compared to that of MgB2. 

The value of A in MgB2, averaged over four groups who 
have done extensive calculations, ||lC|-|l^ is 0.82±0.08. 
The enhancement of 2.1 (above) gives A ~ 1.75 for 
Lio.sBC. We obtain Tc from the Allen-Dynes equation 
with a single frequency; choosing /i* — 0.09 with 
our calculated lu — 58 meV (similar to |l0|-[l3||) and 
A — 0.82 gives Tc = 39 K for the reference material 
MgB2, as observed. Using the same fi* , our calculated 
value iOE2g — 68 meV, and A — 1.75 gives Tc ~ 115 K 
for 2:=0.5. Variation with x should be dominated by the 
a DOS NCT(Ef-(a::)), and in Fig. 5 we plot the calculated 
Tc(a;) assuming other quantities are constant to give an 
idea of the rapidity of the variation of Tc{x). This esti- 
mate indicates that Tc should be higher than 40 K for 
X > 0.25. This estimate neglects multiband effects which 
are substantial for MgB2, [p^ but our purpose of illus- 
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trating the likely value of Tc relative to MgB2 should be 
reasonable. Even if our method of estimating Tc overop- 
timistic, Li3;BC should still be an impressive supercon- 
ductor. 
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FIG. 5. Calculated Tc versus hole concentration, but ac- 
counting only for the x dependence of Nct(Ef). The curve is 
based on the model described in the text that is pegged to re- 
produce MgB2 exactly with A=0.82, ^i* = 0.09, a; = 58 meV. 
The inset shows the A dependence of Tc for MgB2 (calculated 
w = 58 meV) and Lio.sBC (calculated lj = 68 meV). 

According to Worle et al. the hole doping we have 
treated has already been achieved, but measurements at 
reduced temperature were not made. In the event that 
doping is difficult, LiBC (being a semiconductor) seems 
to be a prime candidate for field-effect doping (FED) 
to achieve superconductivity. FED has already been re- 
ported to produce (f ) Tc = 52 K in field-effect hole-doped 
Cgo and at Tc = 117 K in expanded Cgo (2^, and (2) at 
f 4 K in a ladder cuprate and 89 K in an "infinite layer" 
cuprate |^^. Not only do these studies establish the ca- 
pabilities of the FED technique, they demonstrate that 
phonon-coupled superconductivity (as in the fuUerenes) 
can produce very high values of Tc, as we are suggesting 
here for Li^^BC. 

EP coupling of the strength that we obtain for Lia;BC 
must raise the question of the associated lattice insta- 
bilty that will ultimately occur. In MgB2 the i?2g mode 
remains stable although it acquires a very large width 
PI and no doubt is renormalized downward strongly p^ ] 
by the large coupling to the cr states. [In weakly cou- 
pled AIB2 LUe2g = 123 meV. 0] With its much stronger 
coupling, Lia;BC provides the opportunity to approach 
- and perhaps encounter - the associated instabilty. We 
are currently calculating ujE2g i^) to study this instability; 
note however that supposing a (say) 20% lower frequency 
actually leads to a higher Tc in this regime, because the 
larger value of A more than compensates the decrease in 
energy/temperature scale. Raman and inelastic neutron 
spectroscopy will be very informative in revealing the re- 
sponse of the lattice to this large EP coupling. 



To gather these results together: we have shown that 
in comparison to MgB2, the material Lia;BC, x ^ 0.5, 
has (i) the hole-doped a states at E^? that are crucial in 
MgB2, (ii) N(Ef) is comparable to that of MgB2, (in) 
an E2g deformation potential of the a bands that is over 
40% larger than in MgB2, and (iv) the E2g mode fre- 
quency is similar to that in MgB2. In combination, these 
results strongly suggest a coupling strength perhaps as 
high as A ~ 1.75, and a superconducting critical temper- 
ature that could be more than twice as high as in MgB2 
for X = 0.5 ±0.2.. 
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Deutscher Akademischer Austauschdienst. 
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